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Abstract 

q ■ We analyze the impact-parameter dependent diagonal non- singlet generalized 

CN . parton distribution of the pion in two distinct chiral quark models: the Spectral 

+j \ Quark Model and the Nambu-Jona-Lasinio model. Leading-order perturbative 

^ ■ QCD evolution from the low quark-model scale to higher scales is carried out 

with the help of the inverse Mellin transform. The model agrees very rea- 
sonably with the recent results from transverse-lattice calculations, while the 
forward parton distribution agrees with the experimental data. 

> 

00 . 

T^j- . This talk is based on our recent work [ 1 1, where more details and references can be found. Re- 

cently, transverse-lattice calculations have produced first data Q on the impact-parameter dependent 
diagonal non-singlet generalized parton distributions (GPD) of the pion (cf. contribution by Dalley). 
These GPD's have been a subject of intense studies (for a review and literature see, e.g., Ref. 01|4|). The 
impact-parameter (b) formulation, discussed here by Burkardt, is a basis for the transverse lattices 0, 
as shown by van de Sande. In this talk we present predictions for the GPD of the pion from two chiral 
quark models: the recently-proposed Spectral Quark Model (SQM) |6] and the Nambu-Jona-Lasinio 
Qh! model with the Pauli-Villars regulator (NJL). Chiral model results for the skewed 6-integrated (forward) 

GPD were shown here by Vento. 

The off-forward (A_l ^ 0) diagonal (£ = 0) GPD of the pion (we take 7r + ) is defined as 



H(x,Z = 0,-Al) = J d 2 b J ^ e <(^-+Ax-b) {7r + (?/) |^ (0) _L_ )b)7 + 9(0) l_ jb) | 7r + (p))j (1) 

where x is the Bjorken x, and = p' — p lies in the transverse plane. In chiral quark models the 
evaluation of H at the leading- N c (one-loop) level amounts to the calculation of the diagram of Fig. [Ha), 
where the solid line denotes the propagator of the quark of mass u. The calculation is carried in the 
light-cone coordinates and in the Breit frame, A + = ITU , yielding in the chiral limit 

o O 1 I K ± -A ± (l-x) 



7rf4 J (2ny (K_l + (1 - x)A ± y + u 2 

where = (1 — x)p± — xk±. To proceed further, we need to specify the regularization. First, 
we consider the recently proposed Spectral Quark Model (SQM) |6|, described in ERA's talk. The 
approach is very successful in describing both the low- and high-energy phenomenology of the pion, 
while complying to the chiral symmetry and anomalies. In the meson dominance version of the model 
one obtains (D 

9 x m 2 n (m 2 n - (1 - x) 2 A^ ) 
[m A + (1 - x) z Aj_) z 
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Fig. 1: (a) The diagram for the evaluation of the GPD of the pion in chiral quark models, (b) Valence impact-parameter 
dependent diagonal GPD of the pion, V(x,b), plotted as a function of the Bjorken x variable. The results of the chiral quark 
models at the model scale of Q = Qo = 313 MeV. Solid lines: SQM of Ref. 1 6 1, dashed lines: the NJL model with two Pauli- 
Villars subtractions. Label all denotes the forward distribution, i.e., the function V(x,b) integrated over the whole 6-plane. 
Labels [i,j] denote the function V(x,b) integrated over the square plaquettes centered at coordinates (ibo,jbo) of the edge 
of length bo, times the degeneracy of the plaquette (see the text for details). Following Ref. |2|, the value of bo is taken to be 
2/3 fm. (c) The results for V(x,b) at the scale Q ~ 500 MeV, obtained from transverse-lattice calculation of Ref. 0. The 
model results of (b) can be compared to the data of (c) only after a suitable QCD evolution. 

We check that H(x, 0, 0) = 1 and Jq 1 dxH(x, 0, t) = m 2 p / {m 2 p + t), which is the built-in vector-meson 
dominance principle. The NJL model with the Pauli-Villars regularization also yields simple expressions. 

Our aim is to compare our results, after a suitable QCD evolution, to the transverse-lattice data of 
Ref. . These data give the GPD of the pion at discrete values of the impact parameter b, corresponding 
to a square lattice with spacing of bo ~ 2/3 fm. In order to mimic the lattice in our study, the model 
predictions are smeared over square plaquettes, the same ones as in the discrete lattice. The plaquettes 
are labeled [i, j], which means that they are centered at coordinates (ibo,jbo), and have the edge of length 
bo = 2/3 fm Q. The smeared valence GPD is defined as 

V(x, M) = / dh db 2 V( Xl Jb\ + b 2 2 ). (4) 

J(i-l/2)bo J(j-l/2)b V 

Figure [flb) shows the model results with smearing. The degeneracy factor equal to the number of 
plaquettes equidistant from the origin is included, i.e., the [1, 0], [1, 1], and [2, 0] plaquettes are multiplied 
by a factor of four, while [2, 1] would be multiplied by eight. 

Figure 2(c) shows the data from the transverse-lattice calculations shown by Dalley Q. They 
refer to the scale Q ~ 500 MeV. Since the scale pertaining to our calculation is much lower, we need 
to evolve our results upward before comparing to the data of Fig. 2(c). Our model calculation has 
produced distributions corresponding to a low quark model scale, Qo- A way to estimate this scale is 
to run the QCD evolution upward from various scales Qo up to a scale Q where the data can be used. 
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Fig. 2: Results of the LO DGLAP evolution of the impact-parameter dependent diagonal non-singlet GPD of the pion, 
V(x, b, [i,j]), started from the initial condition at Q = Qo = 313 MeV produced by the Spectral Quark Model (Fig. 2(a), 
solid lines). Figures (a,b,c) correspond to Q = 400 MeV, 500 MeV, and 2 GeV, respectively. Labels as in Fig. 2. Figure (d) 
shows xV(x, 6, [i, j]) for Q = 2 GeV, with the dashed line showing the SMRS 1 10 1 parameterization of the data for the forward 
parton distribution function. 



Alternatively, one may use the momentum fraction carried by the quarks at the scale Q and the downward 
QCD evolution in order to estimate Qo- The resulting value is Qo = 313+ioMeV. 

We apply the QCD evolution to the smeared functions of Eq. ©. Figure 2 (a) shows the plaquette- 
averaged functions V(x, Qo, [hj]) f° r SQM (solid lines) and the NJL model (dashed lines). The results 
of the evolution for SQM (for NJL the effect is similar) are shown in Fig. 3 at three values of the 
reference scale Q: 400 MeV (a), 500 MeV (b), and 2 GeV (c). We note a large effect of the evo- 
lution on the distribution functions. The lines labeled all correspond to the forward case, i.e., show 
/ d 2 bV(x, Q, b) = V(x, Q, = 0). The originally flat distribution of Fig. 2(a) recovers the correct 
end-point behavior at x — > 1 according to Eq. ®. As Q increases, the distribution is pushed towards 
lower values of x, as is well known for the DGLAP evolution. At Q = 2 GeV the result agrees very 
well (9) with the SMRS parameterization of the pion structure function ITOl . as can be seen from Fig. 3(d) 
(here we plot for convenience xV(x, Q)) by comparing the dashed and solid lines. The qualitative agree- 
ment of Figs. 1(c) and 2(b) is striking, baring in mind the simplicity of our approach. 

An interesting feature of LO DGLAP evolution is the induced suppression at x —> 1: a function 
which originally behaves as V(x, Qo,b) — > C(b)(l — x) p evolves into ff2l 

_4£p q(Q) 

V{x,Q,b) ^C{b){l-xf log -(Q ), x^l. (5) 

In our approach the integrated function at Q = Qo has p = and the exponent becomes 1.1 ± 0.1 
at Q = 2 GeV and 1.3±0.1atQ = 4 GeV. Note that the Brodsky-Lepage counting rules for the 
behavior at x — > 1 are disobeyed, as well as the predictions are different from the Dyson-Schwinger 
model presented here by Roberts, where the behavior is (1 — x) 2 . On the other hand, our predictions 
agree within experimental uncertainties with the experimental data. Fig. 3 confronts our results evolved 
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Fig. 3: Model prediction for the valence forward parton distribution of the pion, evolved to the scale of 4 GeV, and the E615 
Drell-Yan data US The behavior at x 1 is (1 - x) 1 3±0 \ 

to the scale of 4 GeV to the E615 experimental Drell-Yan data lfT3ll which cover the large-x region. The 
quality of this comparison is impressive. 

In summary, the obtained agreement of our approach, based on non-perturbative chiral quark 
models in conjunction with perturbative LO DGLAP evolution, with the data from the transverse lattices, 
is quite remarkable and encouraging, baring in mind the simplicity of the models and the apparently 
radically different handling of chiral symmetry in both approaches. We also note that the forward parton 
distribution functions are in agreement with the SMRS parameterization and experimental Drell-Yan 
data. 
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